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ABSTRACT 

The age difference between the two main globular cluster sub-populations in the 
Virgo giant elliptical galaxy, NGC 4472 (M 49), has been determined using HST WFPC2 
images in the F555W and F814W filters. Accurate photometry has been obtained for 
several hundred globular clusters in each of the two main sub-populations, down to 
more than one magnitude below the turn-over of their luminosity functions. This allows 
precise determinations of both the mean colors and the turn-over magnitudes of the 
two main sub-populations. By comparing the data with various population synthesis 
models, the age-metallicity pairs that fit both the observed colors and magnitudes have 
been identified. The metal-poor and the metal-rich globular clusters are found to be 
coeval within the errors (~ 3 Gyr). If one accepts the validity of our assumptions, 
these errors are dominated by model uncertainties. A systematic error of up to 4 Gyr 
could affect this result if the blue and the red clusters have significantly different mass 
distributions. However, that one sub-population is half as old as the other is excluded 
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at the 99% confidence level. The different globular cluster populations are assumed to 
trace the galaxy's major star-formation episodes. Consequently, the vast majority of 
globular clusters, and by implication the majority of stars, in NGC 4472 formed at high 
redshifts but by two distinct mechanisms or in two episodes. 

The distance to NGC 4472 is determined to be 15.8 ± 0.8 Mpc, which is in excellent 
agreement with six of the seven Cepheid distances to Virgo Cluster spiral galaxies. This 
implies that the spiral and elliptical galaxies in the main body of Virgo are at the same 
distance. 

Subject headings: galaxies: individual: NGC 4472, galaxies: star clusters, galaxies: 
formation, globular clusters: general 
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1. Introduction 

1.1. Star formation in early type galaxies traced 
by globular cluster formation 

Arguably the two most important open questions 
about the formation of early-type galaxies are when 
they formed their first stars and when they assembled 
dynamically. While the latter question is probably 
best answered by observations at intermediate and 
high redshift, the answer to the former is likely to be 
found from observations in the more local universe. 

A number of studies have estimated the epoch of 
star formation by comparison of observations of the 
diffuse stellar component in nearby and low-rcdshift 
galaxies to stellar population models (e.g. Renzini 
1999 and Bender 1997 for reviews). However, the use 
of diffuse star light has two main disadvantages. First, 
it is very difficult to disentangle different stellar pop- 
ulations and, second, a relatively recent but (in mass) 
unimportant star formation episode can dominate the 
luminosity weighted quantities used to derive the star 
formation history. 

Both these problems are by-passed when globular 
clusters are used to study the main epochs of star for- 
mation. In nearby galaxies, individual globular clus- 
ters appear as point source, and can be character- 
ized by a single age and a single metallicity. Distinct 
populations can easily be identified (e.g. in a color 
distribution of the globular cluster system) and the 
mean properties of the sub-populations can be deter- 
mined. The relative numbers in each sub-population 
will also indicate the relative importance of each star 
formation episode. 

There is strong support for the assumption that 
globular cluster formation traces star formation. First, 
major star formation observed today in interacting 
galaxies is accompanied by the formation of massive 
young star clusters (e.g. Schweizer 1997, although 
their exact nature is still under debate, Brodie et 
al. 1998). Second, the specific frequency (i.e. num- 
ber of globular clusters normalized to the galaxy star 
light, see Harris & van den Bergh 1981) is constant to 
within a factor of two in almost all galaxies (see com- 
pilation in Ashman & Zcpf 1998), which points to a 
link between star and cluster formation. Moreover, 
although this relation seems to fail in giant ellipti- 
cals, McLaughlin (1999) has shown that the number 
of globular clusters is in fact proportional to the to- 
tal mass, and the number per unit mass is extremely 



constant from galaxy to galaxy. Finally, in less vio- 
lent star formation episodes, Larsen & Richtler (1999) 
have shown that the number of young star clusters in 
spirals is directly proportional to the current star for- 
mation rate. 

To first order, globular clusters are likely to be ex- 
cellent tracers for the star formation episodes in early- 
type galaxies. 

1.2. Several epochs of star formation 

There is now evidence for multiple globular clus- 
ter sub-populations in a number of luminous, early- 
type galaxies (Zcpf & Ashman 1993, Ashman & Zepf 
1998 for a recent compilation). The origin of these 
sub-populations is still under debate (e.g. Gebhardt 
& Kissler-Patig 1999) but it seems clear that the main 
sub-populations formed in different star-formation 
episodes, perhaps at different epochs. In the best 
studied galaxies with multiple populations, the num- 
ber of clusters in each sub-population is roughly 
equal, suggesting that these galaxies have had more 
than one major star-formation episode. The over- 
all evidence from photometry and more recently from 
spectroscopy (Kissler-Patig et al. 1998a, Cohen et 
al. 1998) suggests that in the cluster giant ellipti- 
cals the main sub-populations are old. In other non- 
interacting early-type galaxies, the non-detection of 
sub-populations, cither due to their absence or to a 
conspiracy of age and metallicity, leads to the same 
conclusion: no major globular cluster (i.e. star) for- 
mation has occurred since z ~ 1 (Kissler-Patig et 
al. 1998b). 

There have been a few previous attempts to de- 
termine the relative age difference between globular 
cluster sub-populations. In the SO galaxy NGC 1380, 
Kissler-Patig et al. (1997) derived an age difference of 
around 3 to 4 Gyr between the halo and bulge globu- 
lar clusters. In M 87, Kundu et al. (1999) estimated 
the metal-rich clusters to be 3 to 6 Gyr younger than 
the metal-poor ones. As a comparison, in the Milky 
Way, the difference between halo and bulge/disk glob- 
ular clusters seems to be small (less than 1 Gyr: Or- 
tolaniet al. 1995, around 17%: Rosenberg et al. 1999). 

Deriving relative ages from photometry alone is 
complicated by the fact that broad band colors suf- 
fer from an age-metallicity degeneracy (e.g. Worthcy 
1994). One possible solution to this problem is to 
measure several quantities that are affected differently 
by age and metallicity and to combine the results. For 
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example, most broad band colors are more affected by 
metallicity than by age. On the other hand, optical 
magnitudes for a given mass are more affected by age 
than they are by metallicity. A combination of both 
can break the age-metallicity degeneracy, when the 
mass distributions are known (see Sect. 4.1). 

Both a mean color as well as a "mean magni- 
tude" can be measured for a globular cluster sub- 
population. The mean color is simply the mean color 
of the globular cluster sub-sample, and can be de- 
termined from a color distribution of the globular 
clusters. The turn-over (TO) point of the globular 
cluster luminosity function (GCLF) of a given sub- 
population can be used as its "mean magnitude" . The 
absolute TO magnitude of the GCLF has been found 
to be fairly constant from one globular cluster sys- 
tem to the next in a wide variety of galaxy types and 
environments with a peak at M™ ~ — 7.5 (see Whit- 
more 1997 and Ferrarese et al. 2000 for summaries). 
Three quantities affect the TO value: age, metallicity 
and mass. The effects of age and metallicity on the 
magnitude of a single stellar population were quanti- 
fied by Ashman, Conti & Zepf (1995) using popula- 
tion synthesis models from Bruzual & Chariot (1995). 
The mean mass influences the TO magnitude because 
the TO of the GCLF corresponds to a break in the 
globular cluster mass distribution (e.g. McLaughlin 
& Pudritz 1996). That is, if the characteristic mass 
(at which the slope of the power law mass function 
changes) varies, the TO magnitude will be affected. 
However, the mass distributions seem to be extremely 
similar from galaxy to galaxy and are expected to be 
very similar for different sub-populations within the 
same galaxy if destruction processes are the dominant 
mechanism for shaping the mass distribution, and had 
time enough to act on both sub-po pula tions . We will 
come back to this point in Section |4.2| . 

Neither the mean color nor the TO magnitude 
alone can be used to derive the age and metallicity of 
a globular cluster population because of the above- 
mentioned degeneracy. This paper marks the first 
time that both the mean color and the TO magnitude 
of each globular cluster sub-population have been si- 
multaneously compared to s ingle stellar population 
(SSP) models (see Sect. 4T) in order to distinguish 
age and metallicity. Deep high-quality photometry 
is required to reliably determine the mean color and 
the TO magnitude (implying good photometry be- 
low this TO point) and clean globular cluster samples 
are needed, i.e. good discrimination between globu- 



lar clusters and fore-/background contamination is 
essential. The Wide Field and Planetary Camera 
2 (WFPC2) on-board the Hubble Space Telescope 
(HST) provides the most accurate photometric data 
available and makes such an analysis feasible. 

1.3. Globular clusters in NGC 4472 

For this first study, we selected a galaxy that hosts 
a populous globular cluster system which provides 
good number statistics even after the sample is split 
into the different sub-populations. 

NGC 4472 (M 49), the brightest giant elliptical 
galaxy in the Virgo cluster, is located in the center 
of a sub-group within the cluster and is of Hubble 
type E2. It is one of the brightest early-type galaxy 
within 20 Mpc with an absolute magnitude of My ~ 
—22.6 (See Table 1 for other characteristics). It hosts 
6300± 1900 globular clusters (Geisler et al. 1996). The 
specific frequency of globular clusters is Sn = 5.6 ± 
1.7, close to the average for bright cluster ellipticals. 

Previous photometric studies of the globular clus- 
ter system of this galaxy include those of Cohen 
(Gunn-Thuan photometry, 1988), Couture et al. (John- 
son photometry, 1991), Ajhar et al. (Johnson pho- 
tometry, 1994), and Geisler et al. (1996) and Lee et 
al. (1998, both Washington photometry). 

The two most recent studies were the first to ob- 
tain a comprehensive enough data set to allow the 
identification and characterization of the two main 
globular cluster sub-populations. The mean metal- 
licities of the two sub-populations were found to be 
[Fe/H]~ —1.3 and —0.1 dex for the blue and red 
peak, respectively. The observed metallicity gradi- 
ent (A[Fe/H]/A logr ~ —0.4 dex/log(arcsec)) in the 
system is mostly due to the radially varying ratio of 
these two populations. The metal-rich component is 
spatially more concentrated and has the same ellip- 
ticity as the galaxy star light while the metal-poor 
component is more extended and spherical. These 
two components can be loosely associated with the 
"halo" and "bulge" of NGC 4472 suggesting different 
formation mechanisms for these galaxy structures and 
the associated clusters. The relative ages of these two 
components is the subject of this paper. 

The paper is structured as follows: In Section || 
we present the observation and data reduction. Sec- 
tion contains details of the data analysis; color and 
turn-over magnitude determinations, etc. Section 
H compares the results to stellar evolutionary mod- 
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els to determine the age difference between the sub- 
populations. In section ^| we derive a new distance to 
NGC 4472 using the TO magnitude method. Section 
U summarizes the results. 

2. Observations and Data Reduction 

2.1. Observations 

The WFPC2 images of NGC 4472 were obtained 
from two data sets. The nuclear pointing was taken 
from Westphal et al. (from HST program GO. 5236) 
and was complemented by our own pointings to the 
north and south of the galaxy (HST program GO. 5920, 
PI: Brodie) which were offset in DEC by 2' 16" and 2' 
49" from the center of the galaxy, respectively. The 
data are summarized in Tab. 2. 

2.2. Data Reduction 

2.2.1. Image processing 

We used calibrated science images returned by the 
STSclf]. The basic reduction was carried out in IRAF^. 
The images were combined with the task crrej of 
the STSDAS package. The Brodie data set consists of 
paired, dithered exposures, and was thus registered 
by shifting the second of each image pair by 0.506" 
(i.e. 5 pixels for the WF chips, 12 pixels for the PC 
chips) with the task imshift. 

Using the standard aperture radius of Holtzman et 
al. (1995) of 0.5", we tested for any magnitude error 
that might be introduced by incomplete shifting of 
frames. This error was found to be less than 0.0003 
mag for both filters and is therefore negligible. 

2.2.2. Photometry 

The photometry was carried out using the source 
extraction software SExtractor v2.0.19 by Bertin 
& Arnouts (1996). For the purpose of object selec- 
tion only, the images were convolved with the instru- 
mental PSF to improve the detection efficiency. The 

3 Based on observations with the NASA/ESA Hubble Space Tele- 
scope, obtained from the data Archive at the Space Telescope 
Science Institute, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc. under NASA contract 
No. NAS5-26555. 

4 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 



PSF for each filter was taken from the HST Hand- 
book version 4.0, using the PSFs (pixel centered) at 
600 nm for F555W and at 800 nm for F814W. The 
selection criteria were set to 2 connected pixels, 2a 
above the background computed in a 32 x 32 pix 2 
grid. Magnitudes were measured in 2 pixels and 0.5" 
radius apertures (see also next Section). Objects in 
common on the V and / frames were identified and 
matched. The positions on the sky of all objects were 
computed using the task metric in the STSDAS pack- 
age. 

The calibration and transformation to the Johnson- 
Cousins V and / filters followed Holtzman et al. (1995) 
and corrections for charge transfer efficiency and gain 
(7e~ /ADU for all our data) were applied. In what fol- 
lows, the magnitude errors will not include the < 4% 
systematic error potentially introduced by the cali- 
bration. The values measured in the 0.5" radius aper- 
tures were used to compute the final magnitudes in 
order to avoid having to apply an additional correc- 
tion to the data. 

2.2.3. Extended objects calibration 

Globular clusters at the distance of the Virgo 
galaxy cluster are expected to be resolved. Milky 
Way globular clusters have typical half-light radii be- 
tween 3 and 10 pc (e.g. Harris 1996), corresponding 
to 0.04" to 0.13" at 16 Mpc distance. Indeed, most 
of our globular cluster candidates have FWHM mea- 
surements (as returned by SExtractor) slightly larger 
than measured for unresolved objects. This, in turn, 
implies that the photometric corrections proposed by 
Holtzman et al. (1995) for unresolved objects will not 
be valid for converting from the 0.5" radius aperture 
to the total globular cluster magnitude. 

New rough corrections were derived as follows. We 
computed 30 x 30 pix 2 normally sampled PSFs in V 
and / for a WF chip using the Tiny Tim v4.4 pro- 
gram by Krist & Hook (1997), and convolved these 
PSFs with Hubble laws of different core radii (0.05", 
0.1", 0.2", 0.6"). Note that these PSFs include the 
convolution with the diffusion kernel (cf. Krist & 
Hook 1997) which is not well defined. This kernel 
introduces an intrinsic smoothing to the optical PSF 
(as seen at the top of the CCD chip) due to charge 
diffusion. 

The cumulative light distributions of these mod- 
eled objects were compared with those of a set of 
isolated globular cluster candidates of various lumi- 
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nosities and chip positions taken from our data. The 
data in / match almost perfectly the PSF convolved 
with a Hubble law of 0.1" core (corresponding to ~ 8 
pc at a distance of 16 Mpc), while the V data lie in 
between the PSF models convolved with 0.05" and 
0.1" Hubble laws. 

The same procedure, using 10 times over-sampled 
PSFs without applying any diffusion kernel, yielded 
best fitting Hubblc-core radii of 0.3" for / data and 
0.1" for the V data. Until the behavior of the diffu- 
sion kernel is better understood (wavelength depen- 
dencies, etc.), it is possible that a systematic core- 
radius error of about 0.1" (~ 8 pc) may be present 
in all radius estimation studies of Virgo globular 
clusters using WFPC2 data. Nonetheless, we used 
the normally-sampled and diffusion-kernel convolved 
PSFs for our analysis. 

The modeled light profiles (see Fig. |l|), unlike the 
real data, provide high enough signal-to-noise to de- 
termine the extrapolation from an 0.5" radius to the 
total V and I magnitudes. The derived corrections 
are applied to the 0.5" radius aperture measurements 
in addition to the Holtzman point source correction 
of C = 0.1. These additional corrections are Cy = 
0.050± 0.015 and Ci = 0.080± 0.010 (which results in 
a color correction of Ciy_x) = 0.030 ± 0.018). The er- 
rors reflect the uncertainty in matching the data with 
the models and the accuracy of the extrapolation to 
infinity of our fit to the light profile. The errors do not 
include, for example, the fact that the core and the 
half-light radii of different globular clusters may vary 
significantly. In general, our data were calibrated ap- 
plying the correction for point sources of C = 0.10. 
When additional corrections were included, they are 
specifically mentioned in the text. Note that the cor- 
rections are systematic, i.e. they do not affect rela- 
tive magnitudes. Therefore, the accuracy with which 
we are able to determine age-metallicity differences 
should not be seriously affected. 

2.2-4- Artificial star/cluster experiments 

Artificial star/cluster experiments were carried out 
using the tasks starlist and mkobjects from the 
IRAF package, artdata. 1000 objects were added per 
run in both the V and I images. 400 runs were com- 
puted resulting in 400,000 artificial objects per fil- 
ter and covering the whole parameter space of back- 
ground values. Both objects modeled by a WF chip 
PSF and objects modeled by the same PSF but con- 
volved with Hubble laws of core radius 0.05" (V) 



and 0.1" (I) were added. All objects had a color of 
V — I = 1.1, previous tests with 5000 objects having 
shown no significant variations among the complete- 
ness curves for V — I = 0.9, 1.1 and 1.3. 

The selection completeness was computed for ob- 
jects found on both the V and I images. All de- 
tected objects which passed our selection criteria (see 
Sect. 2.2.5) and all input objects were sorted by their 
background value and subdivided in 10 background- 
value bins for V and /, respectively. A skewed Fermi 
function 

F ( X ) = a J^±± (1) 

v ; 1 + cxp(a: - c) d v ' 

was fit to all completeness histograms for each back- 
ground bin. A smoothing procedure using a low-pass 
filter gives nearly the same results. The 50% com- 
pleteness limits lie around V ~ 25.5 and I ~ 24.5 
for the lowest background value. The 50% complete- 
ness magnitude falls off continuously with a growing 
background (see Fig. ^). 

The artificial-object experiments also allowed us to 
determine the error in the photometry and to test the 
accuracy with which the FWHM of objects can be re- 
covered as a function of magnitude. Typical FWHM 
values were deduced for extended objects which are 
a good match to the FWHM values of our globular 
cluster candidates. 

2.2.5. Globular cluster selection criteria 

Object selection criteria were applied to our ex- 
tracted data to provide a clean sample of globular 
cluster candidates for further analysis. Only objects 
which met the color cut of 0.5 < V — I < 1.8 and 
a magnitude cut of V < 25.5 and I < 25.2 were in- 
cluded so that some of the red unresolved galaxy nu- 
clei were rejected (cf. Fig. |J). In addition, a FWHM 
cut in the range 0.75 pix < FWHM < 2.5 pix was 
applied. We also calculated the magnitude difference 
obtained through apertures of 0.5 pix and 3 pix radii 
which allowed us to distinguish between point-source 
objects and extended objects like globular clusters 
(e.g. Schweizer et al. 1996). A very rough upper cut 
of A(0.5-3) > 1.0 was used. 

The neural network of SExtractor v2 . . 19, i.e. the 
CLASS parameter was used to supplement the above 
selection criteria. This parameter was defined for 
well-sampled images (Bertin & Arnouts 1996). Since 
the WFPC2 provides under-sampled images the re- 
liability of the CLASS parameter is reduced. We 
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applied only the selection criterion, CLASS < 1.0, 
to our data. A test showed that for objects brighter 
than the 50% completeness level the CLASS param- 
eter cut reliably eliminates point sources and cosmic 
rays. 

Our final list of globular cluster candidates includes 
705 objects in V and I. This is 38% of the initially 
extracted sample. An electronic list of the globular- 
cluster candidates is available from the authors. 

2.2.6. Background contamination 

One significant advantage of HST data over ground- 
based observations is high spatial resolution which 
allows high-confidence rejection of extended objects. 
However, background contamination is still expected. 
In order to determine the background fraction in our 
data we ran the same extraction and selection criteria 
for the Hubble Deep Field North (HDF-N, Williams 
et al. 1996) images. Since the HDF-N has galactic 
coordinates of I = 125.89 and b = +54.83 while our 
exposures have I = 286.92 and b = +70.20, we might 
overestimate the contamination due to foreground 
stars using this approach, although foreground star 
contamination is, in any case, expected to be negli- 
gible. We obtained from the HST Archive images in 
F606W and in the F814W filter with exposure times 
of 20200 sec for F814W and 15200 sec for F606W, re- 
spectively. These images are far deeper than our own. 
Both magnitudes for F606W and F814W were trans- 
formed to Johnson V and Cousins / magnitudes. We 
determined contamination histograms with and with- 
out the application of our selection criteria in order to 
establish their reliability. Exponential laws have been 
fit to the histograms, after applying the rejection cri- 
teria, up to the 50% completeness limits. The expo- 
nential laws act as input for our Improved Maximum- 
Likelihood Code (see Sect. [3.2.1 ). 

Note the efficiency of our selection criteria. Nearly 
all background objects are removed from the data up 
to the 50% completeness magnitude. The contamina- 
tion is less than 5% for both the Vfsq6W aud Ifsuw 
filters (see Fig. ||). 

We stress again that all the above corrections were 
applied to the full sample, i.e. differential values be- 
tween sub-samples remain unaffected even if small 
systematic effects have been introduced. 



3. Analysis 

Our goal is to determine the mean colors and turn- 
over magnitudes of the luminosity functions of the two 
main globular cluster sub-populations in both the V 
and the / filters. These colors and magnitudes can 
then be compared to the predictions of population 
synthesis models for given ages and metallicities. 

3.1. Mean colors 

In order to determine the mean colors of the sub- 
populations we selected from the total sample objects 
with small photometric errors. The application of a 
magnitude cut at V = 23.75 resulted in errors A(V — 
I) < 0.05. We further applied an upper magnitude 
cut at V = 20.0 (corresponding to My ~ —11.2 at the 
distance of Virgo) in order to exclude any foreground 
star contamination, and a broad cut in color (0.8 < 
(V — I) < 1.8) to isolate the region where globular 
clusters are expected. The color-magnitude diagram 
of our sample is shown in Fig. ^, together with the box 
(dotted lines) of pre-selected objects (see Sect. |2.2.5| ). 
A (V — I) histogram of those selected objects is shown 
in Fig. |. 

Formally, a key underlying assumption is that the 
color (metallicity) does not significantly vary with 
magnitude (mass) within a sub-population. This ap- 
pears to be the case as far as can be tested with our 
data. 

A KMM test (Ashman, Bird & Zepf 1994) was 
run on our constrained sub-sample. The test de- 
termines the confidence with which bimodality can 
be said to be present, the most likely relative con- 
tribution of each population to the sample, and the 
peaks (i.e. mean colors) of each sub-population. We 
find that the color distribution is bimodal at the 
> 99.99% confidence level. The relative contribu- 



tions of objects to the selected sample (cf. Sect. 2.2.5 ) 
are 282 and 423 in the blue and red populations re- 
spectively. The peaks of the color distribution lie at 
(V - I) = 0.99 ± 0.01 and 1.24 ± 0.01 for the blue 
and red sub-populations, respectively. Formally, the 
peaks can be determined with an accuracy of 0.005 
mag. The error adopted here reflects the range of val- 
ues derived by varying the selection criteria (i.e. mag- 
nitude and color cuts). 

Two fitting modes were used: 1) forcing identical 
dispersions (homoscedastic fit), 2) allowing indepen- 
dent dispersions (heteroscedastic fit). We note that 
substantially better fits, in terms of reduced x 2 , are 
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obtained if the dispersions of the blue and red distri- 
butions are allowed to differ, and that the red popula- 
tion seems to have a color distribution twice as broad 
as the blue one (obiuc — 0.06, er rc d — 0.12). 

We estimated the peak colors with yet another 
method that avoids the assumption of a Gaussian dis- 
tribution for each color peak (intrinsic to the KMM 
code). We computed the density distribution in color 
of the objects with the help of a kernel estimator us- 
ing the Epanechnikov kernel (e.g. Gebhardt & Kissler- 
Patig 1999). We used h opt = 0.043 for the smooth- 
ing parameter (see Silverman 1986 for detailed discus- 
sion), an average of the optimal smoothing parame- 
ters for the blue and red peak (/i opt (blue) = 0.023 
and /i pt(red) = 0.063) estimated from the cx-guess 
of KMM for both modes. The density distribution is 
overplotted in Fig. |[ A Biweight fit to each peak re- 
turns (V-I) = 1.00 ±0.01 and 1.22 ±0.01, in perfect 
agreement with the results from the KMM code. 

To obtain the mean colors of the peaks, the above 
(V — 7) values need to be cor rected for the cali- 
bration offset described in Sect. 2.2.3| (A(V — I) = 
0.030 ± 0.018), and for a Galactic reddening towards 
NGC 4472 of E(_B - V) = 0.0224 (Schlegel et al. 
1998), i.e. E(V-J) = l.3-E(B-V) = 0.029 (following 
Dean et al. 1978). Note that the relative color differ- 
ence, later used to derive the age difference, is only 
affected by these corrections in second order. The fi- 
nal, corrected mean colors of the two populations are 
(V - I) = 0.93 ± 0.02 and (V - I) = 1.18 ± 0.02. 

3.2. Turn over magnitudes 

The TO magnitudes of the luminosity functions for 
the blue and red sub-population were individually de- 
termined in each filter using a Maximum Likelihood 
approach (cf. Seeker 1992). The changing background 
counts over our whole sample motivated us to develop 
a version of Seeker's (1992) Maximum Likelihood Es- 
timator which includes background variations. We 
calculated the completeness as a function of magni- 
tude and background noise as described in Sect. 2.2.4 . 
This improved code, described below, was used to de- 
rive the peak of the magnitude distribution for each 
filter and each sub-population. 

Seeker (1992) showed that the GCLF of the Milky 
Way and M31 was best represented analytically with 
the Student's t 5 -function 
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We adopted the Student's t 5 -function as the most 
appropriate distribution for the GCLF of NGC 4472 
since it matches the tails of the distribution better 
than a Gaussian function. The difference in the result 
using the one or the other function is, however, very 
small. 

3.2.1. Improved Maximum Likelihood 

Usually a dataset of n independent observations, 
like a sample of magnitudes M = (mi, mi, to„), is 
distributed in an unknown way. Assuming a distri- 
bution function <p(M\tt) which describes the distribu- 
tion of data properly and depends on q parameters 
fi = (pi,P2) •■•tPq) (particularly the turn-over magni- 
tude to and the dispersion at) the likelihood function 
l(fl\M) — <f>(M\Q) is considered as a function I of 
for a given dataset M. The likelihood function is 



l(Q\M) = Y[(f>(mi\n). 



(3) 



For convenience, during the calculation the logarith- 
mic likelihood function is used. During the evaluation 
the best set of parameters Cl is found if 



dl(Q) , 

on 



|n=n 



0. 



(4) 



For a detailed treatment of Maximum Likelihood the- 
ory see, for instance, Bevington & Robinson (1992). 

Our improved Maximum Likelihood code, based 
on Seeker's (1992) code, takes into account the influ- 
ence of the background-noise level p on the complete- 
ness function I(m,p), in addition to the magnitude. 
Similarly, the photometric error function e(m, p) now 
also depends on these two variables. The distribution 
function is 



4>{m\m , a t , p) = 

K ■ $(m|m , a u p) + (l-K)- B{m\p) 



(5) 



where m a is the turn-over magnitude of the given 
distribution and a t is its dispersion (note: CGauss = 
1.291 • at 5 ), while p is the background noise for the 
evaluation position on the chip. K — Nqq/ '(Nqc + 
A^hdf) is the mixing parameter, i.e. the fraction of 
globular-cluster candidates in the whole data set in- 
cluding the background objects. Note that, in our 
case, K is close to 1. 

It is necessary to normalize all the functions used, 
which is implicitly accomplished by 



(2) 



$(m\m ,a t ,p) = 

I(m\p) ■ (/ t 5 (m'\m , a t ) ■ e(m\m', a t ) dm') 



(G) 
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and 



B(m\p) = I(m\p) ■ b(m) 



(7) 



where b(m) is the background contamination function 
obtained from the HDF data, e is the photometric 
error function 



e(m\m , p) 



1 



exp 



'\2 



(to — to') 

2TT-a(m',p) ^ \2-a{m',p) 



(8) 

which is assumed to introduce Gaussian errors to the 
distribution function <fr with a magnitude and back- 
ground noise dependent dispersion 



a(m, p) = ^W0A-( m ~a) + p 2 . 10 0.8-(m-fc)_ (9) 

The parameters a and b can be estimated with suffi- 
cient accuracy from a two-dimensional fit to the pho- 
tometric errors obtained from the artificial star exper- 
iments. They act as input parameters to our code. 

The need to account for the variable background, 
i.e. the need for the improved code, is best seen in 
Fig. ||. The selection completeness as a function of 
magnitude is plotted for a few typical values of back- 
ground counts found in our datasets. The 50% com- 
pleteness limits shift to brighter magnitudes with in- 
creasing background, i.e. photon noise. 

Alternatively one could approximate a "mean" 
completeness by distributing the artificial stars with 
the same density profile as the globular clusters, with- 
out including the background as a parameter in the 
Maximum Likelihood analysis. The approach we have 
adopted is the more accurate solution to the problem. 
The package described above is available as an IRAF 
task from the authors. 

3.2.2. Maximum Likelihood estimations 

The Maximum-Likelihood routine of Drenkhahn 
(1999) was used to determine the TO magnitudes. 
Our completeness functions fit down to the 50% com- 
pleteness levels (Vtaintest — 25.5, /faintest — 24.5) were 
used as the input. The data were convolved with the 
errors as a function of magnitude (see Eq. ||). The 
input gl obula r cluster data were selected as described 
in Sect. 2.2.5. From KMM mixture modeling we de- 
termined that (V — I) = 1.1 is the best value for 
separating blue and red clusters. For the analysis 
we split the full dataset into the central (Westphal) 
pointing and the combination of the halo (Brodie) 
pointings, since these two datasets had slightly dif- 
ferent exposure times, i.e. slightly different limiting 



magnitudes/completenesses. The TO magnitudes for 
the individual pointings are given in Table 3, together 
with an average for the whole sample. The globular 
cluster luminosity functions are shown in Fig. ^. The 
values in Table 3 still need to be corrected for Galactic 
extinction (A v = 0.069 and A I = 0.0 40) a nd the ad- 
ditional aperture correction (see Sect. 2.2.3), however, 
these corrections do not affect the relative difference 
between the TO magnitudes. 

We varied the input completeness and noise limits 
to the Maximum Likelihood code, the resulting TO 
magnitudes all lie within < 0.15 mag. 

As a consistency check, we attempted to reproduce 
the mean (uncorrected) colors of the sub-populations 
(V - i=0.99±0.01 and 1.24±0.01) by computing the 
differences between the V and I TO magnitudes for 
a given population. We obtained (V — I)to = 1-14± 
0.11 for the blue population and (V — I)to — 1-23 ± 
0.16 for the red population, where the errors are the 
errors of the V and I TOs taken in quadrature. The 
color for the blue population derived in this way is 
somewhat redder than that derived from the color 
distribution but the values agree to within the errors 
for both populations. 

3.2.3. Monte-Carlo simulations 

The error in the TO magnitude directly affects the 
error in the estimation of the age d iffer ence between 
the two sub-populations (see Sect. 4.1). We carried 



out Monte Carlo simulations to check whether our 
relatively modest sample sizes (which were as low as 
118 to 232 objects when divided into inner and outer 
samples) affect the results significantly and to estab- 
lish the associated error. 

The blue and red samples were simulated in the 
V filter using the number of objects in each sample 
and assuming a Gaussian distribution with the pa- 
rameters derived above. These artificial samples were 
input into the Maximum-Likelihood code. The 68% 
confidence levels (lcr) on the TO magnitudes, from 
10 4 iterations, were 0.10 and 0.13 for the blue and 
red samples, respectively. These are in good agree- 
ment with the la errors returned by the Maximum- 
Likelihood estimation. 

3.3. Comparison with previous observations 

The derived colors and TO magnitudes were com- 
pared to previous results from other groups. 
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3.3.1. The mean colors 

Geisler et al. (1996) derived metallicity peaks for 
the blue and red sub-populations from Washington 
photometry. They found peaks at [Fe/H]= — 1.3±0.38 
dex and —0.1 ± 0.38 dex for the blue and red sub- 
samples, respectively. Our colors, derived in Sect. [Q , 
can also be translated into mctallicities. However, 
large errors are introduced by the particular choice 
of a conversion relation of the (V — I) color into 
metallicity. Using the relation given by Kissler-Patig 
et al. (1998a) ([Fe/H]= -(4.50 ± 0.30) + (3.27 ± 
0.32) • (V — I)), which is tuned to reflect the non- 
linear behavior at high metallicity we obtain metal- 
licities of [Fe/H]= -1.45 ±0.43 dex and -0.64 ±0.49 
dex. These values are somewhat lower than, but 
within the errors of, the values derived by Geisler 
et al. (1996). Using the transformation of Kundu & 
Whitmore (1998), [Fe/H] = -5.89 + 4.72 • (V- 1), our 
data yield metallicities of [Fe/H] = —1.50 ± 0.05 dex 
and [Fe/H] = — 0.32±0.05 dex. The error corresponds 
to our photometric error in (V — I) only. No calibra- 
tion error of Kundu & Whitmore is considered. The 
metallicity calibration of Couture et al. (1990), [Fe/H] 
= —6.096 + 5.05 • (V — I), gives the somewhat larger 
values of [Fe/H] = -1.39 ±0.05 dex and -0.11 ±0.05 
dex, similar to those derived by Geisler et al. (1996). 

3.3.2. The TO magnitudes 

TO magnitudes have been derived by several groups 
but only for the whole sample (i.e. an average of the 
red and blue populations) and in different filters. Pre- 
vious results should be compared to the average of our 
blue and red sample, with the necessary correction for 
reddening and aperture size: V = 23.76 ± 0.20 (in- 
cluding the error on the TO magnitude and the error 
on the aperture correction). Furthermore, all other 
groups assumed E(B - V) = 0.00 towards NGC 4472 
(using Burstein & Heiles 1982/84) while we adopted 
the Schlegel et al. (1998) map from DIRBE/IRAS 
whose zero-point is offset by 0.020 in E(£> — V). Our 
value would translate into Vq = 23.83 ± 0.20 if we 
assumed E(B -V) = 0.00. 

Using the Washington-Johnson transformation of 
Geisler (1996), V ~ Ti +0.5, Lee et al.'s (1998) result 
translates into Vq = 23.81 ± 0.07, in excellent agree- 
ment with our result. We note that Lee et al. investi- 
gated their blue and red samples separately, but did 
not find any significant difference in their TO magni- 
tudes to within ~ 0.15 mag. However, their analysis 



was less rigorous in this respect than ours, and their 
data suffered from more background contamination, 
as is expected for ground-based observations. 

Harris et al. (1991) derived B T o = 24.78 ± 0.22 
with a 3-parameter Gaussian fit to the GCLF. Using 
B — V ~ 0.8 for the total globular clusters population, 
this translates into Vto = 23.98 ± 0.22. Ajhar et 
al. (1994) derived R TO = 23.3 ± 0.2 (read off their 
Fig. 14) which, for V — R ~ 0.5, translates into Vto = 
23.8 ± 0.2. These results are in very good agreement 
with our determination. 

3.4. Radial Dependencies 

3.4-1- TO magnitudes as a function of radius 

Variation in the TO magnitude with radius within 
a population can occur if the mean age, metallicity 
or mass changes with radius. Large age and metal- 
licity gradients are ruled out by the absence of any 
significant color gradient in the sub-populations (see 
Fig. 0). A systematic change in the mean mass 
with radius could be caused by destruction processes 
shaping the globular cluster mass function. These 
are more important towards the center of the galaxy 
(e.g. recent work by Vesperini 1997 and Gnedin & 
Ostriker 1997), and significant changes are only ex- 
pected within ~ 5 kpc of the center (~ 65" at 16 
Mpc). There the peak of the GCLF is expected to be 
brighter by up to ~ 0.3 magnitudes, and the peak 
should be sharper as a consequence of the prefer- 
ential destruction of low-mass clusters. Our inner 
sample (Westphal data) includes clusters at projected 
radii between 2" and 127" with a mean of 65" and 
our outer sample (Brodie data) includes clusters at 
projected radii between 76" and 252" with a mean 
of 160". Differences are therefore expected between 
these two samples. For the blue population, the TO 
magnitude varies from the inner to the outer regions 
by 0.12±0.20 and -0.29+0.25 in V and / respectively. 
The weighted mean is — 0.04±0.16, i.e. we observed no 
brightening of the TO towards the center, although 
the uncertainty is on the order of the expected varia- 
tion. For the red population the TO magnitude varies 
from the inner to the outer region by 0.33 ± 0.17 and 
0.14±0.25 in V and / respectively. Formally, the TO 
appears brighter by 0.27 ± 0.14 (weighted mean) for 
the inner red population. We do not see any changes 
in the dispersion of the GCLF to within the errors. 

This result is probably not yet secure enough to 
be worth commenting on in depth. However, a dif- 
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ference between the red and blue population could 
occur, for example, if the blue and red clusters are on 
different orbits. Sharpies et al. (1998) found tentative 
evidence for the blue clusters in NGC 4472 having a 
higher velocity dispersion and higher rotation than 
the red ones. Kissler-Patig & Gebhardt (1998) found 
a similar (but firmer) result for M 87, another giant 
elliptical that hosts enough globular clusters for such 
an analysis. If indeed the blue clusters are preferen- 
tially on tangential orbits, while the red clusters are 
preferentially on radial orbits, destruction processes 
would be more efficient for the red population. 

Finally, we note that Harris et al. (1998) and 
Kundu et al. (1999) also looked for this effect in M 87. 
Kundu at al. even looked for trends in the blue and 
red populations separately. Neither group found any 
brighting of the GCLF towards the center of M 87, 
but note again that the uncertainties in the measure- 
ments are still of the order of the expected signal. 

3.4-2. Color as a function of radius 

We looked for color gradients within the sub- 
populations using a weighted least squares fitting rou- 
tine. In addition to the color cuts (0.5 < V — I < 1.1 
and 1.1 < V — I < 1.8), the samples were cut at V 
= 24.8 to provide at least 50% completeness even in 
the central regions. The results are plotted in Fig. 0, 
together with a fit for the total sample. The de- 
rived gradients are A(V - I)/ A logr = 0.015 ± 0.004 
mag/log(arcsec) and 0.013±0.007mag/log(arcsec) for 
the blue and red population respectively. Different 
color and magnitude cuts scatter these results within 
0.01 mag/log(arcsec). Any positive gradients in the 
sub-populations are either very small or non-existent. 
For the whole sample, however, we find a clear neg- 
ative gradient of A(V - J)/Alogr = -0.029 ± 0.07 
mag/log(arcsec), corresponding to a metallicity gradi- 
ent of A [Fe/H] / A log r ~ -0.15±0.02 dex/log(arcsec) 
(somewhat dependent on the color-mctallicity rela- 
tion used). 

We confirm Geisler et al.'s (1996) result, that this 
gradient is caused by the changing ratio of red to blue 
clusters as a function of radius, rather than an over- 
all decrease in metallicity with radius. Quantitatively, 
we find a somewhat smaller gradient than reported by 
Geisler et al. (1996), who found A [Fe/H] /A log r ~ 
—0.4 dex/log(arcsec). This is probably due to our 
smaller radial coverage and our use of a less metallicity- 
sensitive color. Indeed, our sample extends to only 
250" while Geisler et al.'s data extend out to ~ 500". 



The ratio of blue to red clusters varies by only ~ 10% 
from our inner to outer regions, as derived from a 
comparison of their cumulative radial distributions. 

A similar situation is seen in M 87, for which Har- 
ris et al. (1998) report no color gradient within 60", 
and Kundu et al. (1999) report a weak gradient of 
A(V - I) /A logr = -0.017 ± 0.012 mag/log(arcsec) 
within 100". However, Lee & Geisler (1993) report 
a strong gradient (A[Fe/H]/A log r = -0.65 ± 0.17 
dex/log(arcsec)) over 500" radius. 

3.5. Globular Cluster Sizes 

The high resolution of HST allows us to study the 
relative sizes (see Sect. 2.2.3] ) of the globular clus- 
ter candidates using the parameter A(0.5 — 3) (see 
Sect. 2.2.5). For extended objects the magnitude dif- 
ference in two different apertures will appear large as 
the light profile contributes a non-negligible amount 
of light to the outer aperture. The opposite is the case 
if the light profile is narrow. For this analysis, we re- 
stricted ourselves to clusters found on the Wide Field 
(WF) camera chips, to avoid the problematic con- 
version between the WF and the Planetary Camera 
(PC). Fig. H presents the histogram of globular cluster 
sizes for both the blue and red sub-populations. The 
median values for these two distributions point to a 
systematic size difference between the blue and the 
red globular clusters. The median values for the blue 
and red populations are toa = 2.06 ± 0.01 mag and 
toa = 2.02 ± 0.01 mag, respectively, indicating that 
the blue globular clusters are substantially larger than 
the red globular clusters. A Maximum-Likelihood es- 
timation of the mean and its error yields the same 
results. 

Interestingly, in both NGC 3115 (Kundu & Whit- 
more 1998) and M 87 (Kundu et al. 1999), the blue 
clusters also appear significantly larger than the red 
ones. The authors suggest this is due to different ra- 
dial distributions for the red and blue clusters, noting 
that Galactic globular clusters with larger galacto- 
centric distances have larger half-light radii (van den 
Bergh 1994). 

To study the radial dependence of globular clus- 
ter size, we plotted the size parameter (A(0.5 — 3)) 
versus radius and found no gradients: A(0.5 — 3) = 
2.07(±0.03) - 0.5(±3.2) • 10~ 4 • r where r is the radius 
in arcsec, and A(0.5 - 3) = 1.93(±0.03) - 2.9(±2.7) • 
10 -4 • r, for the blue and red clusters respectively (see 
also Fig. H). The difference between blue and red 
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clusters is present at all radii. We therefore speculate 
that the different sizes are a relic of the different for- 
mation processes, unless the red and blue clusters are 
on significantly different orbits. For example, the red 
clusters could be preferentially on radial orbits and be 
systematically more affected than the blue clusters by 
the central potential during their close passage near 
the galaxy center. 

4. The age difference between the main sub 
populations 

4.1. Comparison with the models 

We compare our results (TO magnitude and mean 
color) to various population synthesis models in order 
to derive the age difference between the blue and red 
globular cluster sub-populations. 

We used the new models from Maraston (1998) as 
well as the most recent models from the Gottingen 
Group (Kurth ct al. 1999), and the models from 
Bruzual & Chariot (1996) and Worthey (1994). In 
each case we compared our data to every available 
Initial Mass Function (IMF), i.e. usually the Salpeter 
IMF and a multi-slope IMF. The full range of metal- 
licities was used as well as the necessary age range. 
Since we are only interested in age differences the ex- 
act normalization of the absolute magnitude, which 
is dependent on mass and mass-to-light (M/L) ra- 
tio for a given IMF, is not important. The grids in 
Fig. [| to [12] can, therefore, be somewhat arbitrarily 
shifted in the y-direction. In the following, we fixed 
the grids in the y-direction such that the oldest sub- 
population lies on the oldest computed isochrone (15 
Gyr or 16 Gyr). Further, we arbitrarily set the blue 
population to V = and I = in order to compare 
the different models directly. There is no freedom in 
the x-direction since the color is not mass-dependent, 
and the grid is fixed in V — I for a given stellar pop- 
ulation. 

For the mean colors of the blue and red sub- 
populations we used the colors derived in Section [O . 
The TO values used were the ones derived for the full 
sample (see Table ||). We stress again that additive 
corrections to both the blue and red TO magnitudes 
(i.e. reddening correction, potential systematic errors, 
aperture corrections) only influence the difference as 
second order effects, if at all, and are not a concern. 

However, we make two important underlying as- 
sumptions. We assume that the M/L ratios for blue 
and red clusters are only dependent on age and metal- 



licity (i.e. the IMFs of the blue and red clusters are 
similar). Second, we assume that blue and red cluster 
populations have the same underlying mass distribu- 
tion (by which we mean globular cluster masses, not 
IMFs). Under these assumptions, TO magnitudes are 
influenced only by age and metallicity. We briefly dis- 
cuss the validity of these assumptions and the errors 
they potentially introduce. 

Jf.. 1.1. The M/L and mass distribution assumptions 

Kissler-Patig et al. (1998b) review the empirical 
reasons for expecting the mass distributions and M/L 
ratios to be similar among various sub-populations. 
Briefly, the M/L assumption is essentially the as- 
sumption that both blue and red clusters have similar 
IMFs, at least below 1 Mq where it influences our re- 
sult. Observations so far for globular clusters in the 
Local Group seem to confirm this picture (e.g. Du- 
bath & Grillmair 1997). In the case of the mass dis- 
tribution, globular cluster formation seems to follow 
a universal mechanism (e.g. Elmcgrecn & Efremov 
1997). The constancy of the TO magnitude from 
galaxy to galaxy (e.g. the review by Whitmore 1997) 
and the observed mass distributions of young clusters 
(e.g. brief review by Schweizer 1997) support this pic- 
ture. 

The possibility that size differences might indicate 
different formation processes for blue and red clusters 



was discussed in Section 3.5. It is not clear, however 



whether it is only the mean sizes or also the mean 
masses that are affected. Moreover, our assumption 
would fail if, despite similar mass distributions at for- 
mation, blue and red clusters are affected differently 
by destruction processes. There is a hint (see Section 
3.4 ) that this might be the case in the inner region, 
where the difference between the blue and the red TO 
magnitudes is smaller by 0.32±0.24 (combining the V 
and / information) than in the outer region. (Wc note 
for completeness that an age gradient in the popula- 
tion would have the same observable effect as a mass 
gradient). If this effect is real, the TO difference of 
the total sample, as used, could be affected by a 0.1 
to 0.2 mag systematic error. We estimate that this 
empirical determination gives about the right order 
of magnitude for the potential systematic error in our 
method and, in what follows, we will assume a poten- 
tial systematic error of ~ 0.2 mag, which translates 
into roughly 4 Gyr. 

Note that for similar M/L values of the blue and 
red clusters, such a difference of 0.2 mag would trans- 
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late into only a 20% difference in the mass character- 
izing the TO. Irrespective of the accuracy needed for 
our study, blue and red clusters have roughly similar 
mass distributions. 

4-1.2. Avoiding assumptions about the masses 

For completeness, we mention a method for check- 
ing the validity of these assumptions and a method 
for avoiding them. 

For the chosen filters, age and characteristic mass 
influence the TO far more than metallicity, which is 
fairly well determined by the color. Therefore, one 
could determine the mean mass difference (and check 
the above assumption) by obtaining an independent 
estimate of the mean age difference. An independent 
age estimate could be obtained with spectroscopy of 
a large number of globular clusters and the measure- 
ment of an age-sensitive absorption line index. Sev- 
eral tens of H/3 measurements in each sub-population 
would allow the mean H/3 value to be determined to 
within < 0.1 A which would translate into an accu- 
racy in the mean age difference of ~ 2 Gyr. This in 
turn would allow us to verify that the characteristic 
TO masses of the red and blue populations lie within 
10% to 20% of each other. 

An alternative method for determining the mean 
age difference without making any assumptions about 
the mass distributions would be to use mass-indepen- 
dent quantities that depend only on age and metal- 
licity. Broad-band colors are such quantities. A plot 
similar to Figures ^ through |lj, but in a color-color 
plane, would remove any uncertainties due to differ- 
ent mass distributions. However, it is impossible to 
find two colors in the optical and NIR spectrum that 
disentangle age and metallicity as efficiently as the 
color/TO magnitude combination. 

4- 1.3. Comparison with Maraston (1998) 

In Figure ^, we compare our results to the mod- 
els of Maraston (1998). She computed single burst 
population models with a Salpeter (1955) and multi- 
slope IMF (Gould et al. 1998) for metallicities Z= 
0.001,0.006,0.02 and 0.04, and ages up to 15 Gyr. 
We refer the reader to the original paper for further 
details of the modeling. 

The derived age difference, which is quoted in 
the following as (agebiuc)— (age ro( j), varies between 

0.3 Gyr (Gould et al. IMF, I band) and - +2.4 

Gyr (Gould et al. IMF, V band), with errors of about 



2.7 Gyr (la). Within the errors, the blue and red 
population appear coeval, with a formal difference of 
+0.7 ± 1.8 Gyr (straight mean of the four computed 
differences with their dispersion shown as the error). 
A 3<T limit on the age difference would be around 6 
Gyr. 

4-1-4- Comparison with Kurth et al. (1999) 

In Figure |l^, we compare our results to the re- 
cent models of Kurth et al. (1999). They computed 
single burst models with a Salpeter (1955) and a 
Miller-Scalo (1979) IMF for 6 metallicities between 
Z= 0.0001 and 0.05, and ages up to 16 Gyr. 

The age difference varies between —3.8 Gyr (Sal- 
peter IMF, I band) and -0.4 Gyr (Salpeter IMF, V 
band) with a typical error of 2.2 Gyr. The mean 
difference is —2.3 ± 1.6 Gyr, i.e. the red population 
appears somewhat older than the blue one. The 3<r 
differences are around —7 and +2.5 Gyr. 

4-1.5. Comparison with Bruzual & Chariot (1996) 

In Figure [0], we compare our results to the models 
of Bruzual & Chariot (1996). We used the single burst 
populations computed with a Salpeter (1955) IMF 
and a Scalo (1986) IMF, with metallicities between 
Z= 0.0004 and 0.05. Their models with a Miller- 
Scalo (1979) IMF were not used. The tracks for their 
highest metallicity (Z= 0.1) are not shown. The old- 
est modeled populations is 16 Gyr. 

The age differences vary between +1 Gyr (Scalo 
IMF, / band) and +6 Gyr (Salpeter IMF, V band), 
with an typical error of 2.8 Gyr. The mean difference 
is +3.3 + 2.2 Gyr. The 3<r limits lie around —3 and 
+10 Gyr. 

4-1.6. Comparison with Worthey (1994) 



In Figure 12, we compare our results to Worthey's 
(1994) single burst models. He computed models with 
a Salpeter (1955) IMF and a Miller-Scalo (1985) IMF 
for metallicities between Z= 0.0002 and 0.05, up to 
an age of 16 Gyr. 

The age difference varies between —5.8 Gyr (Salpeter 
IMF, / band) and -2.1 Gyr (Miller-Scalo IMF, V 
band), with a typical error of 2.5 Gyr. The mean 
difference is —4.0+ 1.6 Gyr, i.e. the red populations 
appears older than the blue one. The 3cr limits are 
around —9 and +1 Gyr. 
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4- 1.7. The maximum age difference between the glob- 
ular cluster populations 

We note that the I band results give a system- 
atically smaller difference (in the sense we defined 
it: agebiuc— age re d) than the V band results. Except 
for the Bruzual & Chariot (1996) models, the results 
from the / band suggest that the red population is 
marginally older than the blue one. 

The different models predict mean differences any- 
where between —4.0 Gyr and +3.3 Gyr, with a mean 
of —0.6 Gyr and a dispersion of 3.2 Gyr. It would 
be somewhat arbitrary to assign more weight to one 
model than to another, so that we are currently lim- 
ited in the age determination by the model uncertain- 
ties. Nevertheless, it seems that both populations are 
coeval within the errors, and we can exclude at the 
99% confidence level one population being half as old 
as the other. 

If our results are affected by a systematic error, 
e.g. if the smaller sizes of the red globular clusters do 
indeed translate into a lower "mean" mass (such that 
the red TO magnitude is ~ 0.2 lower - see Section 
|4 . 1 . 1| ) . the age difference would be affected by ~ +4 
Gyr. 

4.2. Implications for the star formation episodes 
of NGC 4472 

The above results suggest that the vast majority 
of globular clusters and stars formed at high redshift 
in this cluster galaxy. As already mentioned in Sec- 
tion [I], several other lines of evidence lead to similar 
results for cluster galaxies. Interestingly, the globular 
cluster results imply that two major star-formation 
events took place at high redshift. Those two events 
appear to have happened moderately closely in time 
but from differently enriched gas. The blue clusters 
appear to have formed from gas with less than l/20th 
solar metallicity while the red clusters formed out of 
gas with close to solar metallicity. Furthermore, the 
mean sizes of the blue and red clusters appear to be 
significantly different. This suggests different condi- 
tions at the formation epoch, unless they are directly 
related to the metallicity and associated cooling pro- 
cesses in the molecular clouds in which the clusters 
originated. 

It is unclear whether the blue or the red clusters 
are the older ones. At face value it seems odd that 
the more metal-rich clusters would be the older ones. 
However, in scenarios where the blue clusters formed 



in halo clouds (e.g. Kissler-Patig et al. 1998b) or were 
accreted with dwarf galaxies (e.g. Cote et al. 1998, 
Hilker et al. 1999a, b), this could be understood in 
terms of the large structures (e.g. bulges) collapsing 
first, followed shortly thereafter by the smaller struc- 
tures (clouds/dwarf galaxies). This, however, would 
force a revision of the merger scenario which clearly 
predicts that the red clusters will be younger than the 
blue ones. 

Alternatively, if two sub-populations of clusters 
are coeval or if the blue clusters are the older ones, 
age differences of a few Gyr between the different 
sub-populations could be explained in all globular 
cluster /galaxy formation scenarios, i.e. the above- 
mentioned scenarios, merging (Ashman & Zepf 1992), 
or in-situ formation of distinct populations (e.g. Forbes 
et al. 1997, Harris et al. 1998, Harris et al. 1999). We 
note only that, in the merger scenario, merging would 
have to have taken place at early times (z > 1 at 
the 3cr level) which is somewhat at odds with the re- 
cent theoretical predictions from hierarchical cluster- 
ing models (e.g. Kauffmann 1996, Baugh et al. 1998). 
These tend to predict that the main star formation in 
cluster environments occurred at redshifts z < 1. 

5. A distance estimation for NGC 4472 

The TO magnitude (Vro) of the GCLF appears to 
be essentially universal and, therefore, can be used 
as a distance indicator (e.g. the review by Whit- 
more, 1997). Vro depends on the age and metal- 
licity of the globula r clu sters (Ashman, Conti & Zepf 
1995, see also Sect. 4.1), and weakly on galaxy type 
(Harris 1997) which may be a consequence of the 
age/metallicity dependence. The derived TO mag- 
nitude for the globular clusters in NGC 4472 can be 
used to obtain a distance to the galaxy. 



5.1. The absolute turn over magnitude 

The method is best calibrated using the TO mag- 
nitude of the Milky Way GCLF, since it involves the 
minimum number of steps on the distance ladder. The 
M31 GCLF serves as a strong check on this calibra- 
tion (cf. Ferrarese et al. 2000; Kavelaars et al. 1999). 

The Milky Way system is dominated by old glob- 
ular clusters of mean metallicity [Fe/H]~ —1.4. It 
is, therefore, well-suited to fitting the TO magnitude 
of the blue population in NGC 4472 (see Sect. fi~i"| ) 
and no correction for age-metallicity dependence need 
be applied. However, we have to assume that age- 
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metallicity dominates any dependence on galaxy type, 
and that globular cluster mass functions are intrin- 
sically similar, which seems justified to first order 
(e.g. Kissler-Patig et al. 1998b). 

Delia Valle et al. (1998) recently re-derived Vto 
for the Milky Way system. To derive individual dis- 
tances to all Galactic globular clusters, they adopted 
a horizontal branch (HB) magnitude metallicity re- 
lation derived from globular cluster main sequence 
fitting using distance measurements to sub-dwarfs 
from HIPPARCDS (in this case by Gratton et al. 1997). 
From the resulting GCLF they obtained Vto = 
—7.62 ± 0.06 noting that, depending on the adopted 
HB magnitude-metallicity relation, this result could 
suffer from a to —0.2 mag systematic error. 

Drenkhahn & Richtler (1999) used a different ap- 
proach. They assumed the LMC distance to be the 
fundamental zero-point (here (to — M) = 18.46 ± 
0.06), and derived an RR Lyrae magnitude-metallicity 
relation. The remaining steps were similar to those 
of Delia Valle et al.. Drenkhahn & Richtler found 
V TO = -7.61 ± 0.08 and I TO = -8.48 ± 0.10 for the 
Milky Way GCLF. 

5.2. The distance to NGC 4472 

The apparent TO magnitude derived for the blue 
population is V T o = 23.62 ± 0.09 and I T o = 22.48 ± 
0.07 (see Tab. ||). These values need to be corr ected 



2.2.3 



for the photometric term discussed in Sect 
(ACV = 0.050 ±0.015, AC i = 0.080 ±0.010), as well 
as for Galactic extinction (A v = 0.069, Ai = 0.040). 
Using the absolute TO magnitudes discussed above, 
we derive distance moduli of (m — M)v = 31.11±0.19 
and (to — M)i = 30.84 ± 0.14 (adding all errors in 
quadrature), i.e. a mean (m — M) — 30.99 ± 0.11 
(where the error does not include any possible sys- 
tematic error in the calibration). The corresponding 
distance is 15.8 ± 0.8 Mpc. 

A quick comparison with distances derived for 
NGC 4472 by other methods shows very good agree- 
ment (despite completely different calibrations). Ja- 
coby et al. (1990) derived (m-M) = 30.84±0.11 from 
the planetary nebulae luminosity function. Tonry et 
al. (1990) and Neilsen (1999) derived (to - M) = 
30.78 ± 0.07 and (m - M) = 30.93 ± 0.08, respec- 
tively, from surface brightness fluctuations. 

Cepheid distances now exist for 7 spiral galaxies in 
the Virgo Cluster (cf. Pierce et al. 1994 and Ferrarese 
et al. 2000). With the exception of NGC 4639 (Saha 



et al. 1996), which appears to be background to the 
main body of spirals, the mean distance to the other 
6 galaxies is (to - M) = 31.01 ± 0.07 or 16.0±0.6 
Mpc. Our GCLF distance to NGC 4472 is in excellent 
agreement with this distance, which supports the view 
that the Virgo cluster spirals and ellipticals are at 
essentially the same distance. 

6. Summary 

The age difference between the two major globular 
cluster sub-populations in an early-type galaxy has 
been determined. We have developed an age differen- 
tiating method which uses the precise TO magnitudes 
and colors of the sub-populations. We have applied 
this method to NGC 4472, matching our observables 
to several SSP synthesis models. 

An extended-object analysis showed that at the 
distance of Virgo (i.e. ~16 Mpc) HST is capable of 
resolving single globular clusters. Thus it is necessary 
to make an additive correction (of Cy — 0.050± 0.015 
mag and Cj = 0.080 ± 0.010 mag at the distance of 
Virgo) to the mean correction of C = 0.1 mag (for 
point sources) when extrapolating HST photometry 
obtained in 0.5" radius apertures to total luminosities. 
This corresponds to a color correction of C(y—i) = 
0.03 ±0.018. 

We found the globular cluster color distribution 
in NGC 4472 to be clearly bimodal (at the 99.99% 
confidence level) with two peaks located at V — I = 
0.99±0.01 and 1.24±0.01 corresponding to (V-I) = 
0.93±0.02 and (V-I) = 1.18 ± 0.02 when corrected 
for reddening and finite aperture. We determined 
that the red globular clusters have a color distribu- 
tion almost twice as broad as the blue ones. 

We divided the whole globular cluster population 
into an inner and an outer sample, as defined by 
the HST different pointings. We used an improved 
Maximum-Likelihood estimator, which accounts for 
varying background values, to determine the TO mag- 
nitudes for each sub-population in each radial sample 
(see Tab. ||). The TO of the blue sub-population is 
the same for the inner and outer sample. By contrast, 
the red sub-population appears to be marginally 
brighter by 0.27 ± 0.14 mag in the inner region than 
in the outer region. This difference is of the expected 
order of magnitude if destruction processes act on the 
red sample within 5 kpc. 

The radial color distributions of each sub-population 
show almost no gradients. However, we detected an 
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overall gradient of A(V - J)/Alogr = -0.029 ± 0.07 
mag/log(arcscc) in the system, which is mostly an 
effect of the changing ratio of red to blue globular 
clusters with radius (cf. Geisler et al. 1996). 

Taking advantage of HST's high spatial resolution, 
we estimated the mean sizes of the globular clusters 
using the A(0.5 — 3) parameter. The blue globular 
clusters were found to be larger than the red ones, 
as was already observed in M 87 (Kundu et al. 1998) 
and NGC 3115 (Kundu & Whitmore 1997). The most 
likely explanations are that the size difference reflects 
different formation processes, or that red clusters are 
on preferentially radial orbits, while blue clusters are 
on preferentially tangential orbits. 

We used 4 different SSP models to derive the 
mean age difference between the red and blue sub- 
population. Our underlying assumption was that 
both blue and red populations have similar mass dis- 
tributions which appears to be valid at least at the 
10% to 20% level (which, nevertheless, corresponds 
to a potential systematic error of ~ 4 Gyr). The 
mean age difference between the blue (metal-poor) 
and the red (metal-rich) sub-populations is aget; ue — 
age re d = — 0.6±3.2 Gyr, where the error is the disper- 
sion introduced by the different models (for a given 
model the error on our method is ^2 Gyr). 

We conclude that the globular clusters (and by as- 
sociation the stars) in NGC 4472 formed in two major 
episodes/mechanisms that were coeval within a few 
Gyr. 

The images as well as electronic lists of the photom- 
etry are available from the authors. Also available is 
the maximum-likelihood code as an IRAF task used to 
derive the globular cluster luminosity functions tak- 
ing into account the background variations over the 
field. 

Acknowledgments 

We are very thankful to Claudia Maraston for cus- 
tomizing part of her models for our purposes and pro- 
viding electronic files, as well as to Gcorg Drcnkhahn 
for his help in implementing his Maximum Likelihood 
code. We also thank Sandy Faber, Duncan Forbes, 
Karl Gebhardt, and Klaas de Boer for interesting dis- 
cussions. 

THP gratefully acknowledges support from the 
UCO /LICK Observatory. MKP was partly supported 
by the Alexander von Humboldt Foundation. This 
research was funded by HST grants GO.05920.01- 



94A and GO. 06554. 01-95A and faculty research funds 
from the University of California at Santa Cruz. This 
research was partly supported by NATO Collabora- 
tive Research grant CRG 971552. 

REFERENCES 

Ajhar, E.A., Blakeslee, J.P., Tonry, J.L. 1994, AJ, 
108, 2087 

Ashman, K.M., Zcpf, S.E. 1992, ApJ, 384, 50 

Ashman, K.M., Zepf, S.E. 1998, Globular Cluster 
Systems, Cambridge University Press 

Ashman, K.M., Bird, CM., Zepf, S.E. 1994, AJ, 108, 
2348 

Ashman, K.M., Conti, A., Zepf, S.E. 1995, AJ, 110, 
1164 

Baugh, CM., Cole, S., Frenk, C.S., Lacey, C.G. 1998, 
ApJ, 498, 504 

Bender, R. 1997, ASP Conf. Ser. Vol. 116, The 
nature of elliptical galaxies, eds. M. Arnaboldi, 
G.S. Da Costa, P. Saha 

Bertin, E., Arnouts, S. 1996, A&AS, 117, 393 

Bevington, P.R., Robinson, D.K. 1992, Data Reduc- 
tion and Error Analysis for the Physical Sciences, 
McGraw-Hill, 2nd edition 

Brodic, J. P., Schroder, L.L., Huchra, J. P., Phillips, 
A.C, Kissler-Patig, M., Forbes, DA. 1998, AJ, 
116, 691 

Bruzual, A.C, Chariot, S. 1996, electronically avail- 
able see: Leitherer, C, et al. 1996, PASP, 108, 
996 

Burstein, D., Heiles, C. 1982, AJ, 87, 1165 
Burstein, D., Heiles, C. 1984, ApJS, 54, 33 
Cohen, J. G. 1988, AJ, 95, 682 

Cohen, J.C, Blakeslee, J.P., Ryzhov, A. 1998, ApJ, 
496, 808 

Cote, P., Marzke, R.O., West, M. 1998, ApJ, 501, 554 

Couture, J., Harris, W.E., Allwright, J.W.B. 1990, 
ApJS, 73, 671 

Couture, J., Harris, W.E., Allwright, J.W.B. 1991, 
ApJ, 372, 97 

Dean, J.F., Warren, P.R., Cousins, A.W.J. 1978, MN- 
RAS, 183, 569 

Delia Valle, M., Kissler-Patig, M., Danziger, J., Storm, 
J. 1998, MNRAS, 299, 267 



16 



Drcnkhahn, G. 1999, Diploma Thesis, University of 
Bonn 

Drenkhahn, G., Richtler, T. 1999, A&A, in press 

Dubath, P., Grillmair, C.J. 1997, A&A, 321, 379 

Elmegreen, B.G., Efremov, Y.N. 1997, ApJ, 480, 235 

Ferrarese, L. et al. 2000, ApJ in press 

Forbes, DA., Brodie, J.P., Grillmair, C.J. 1997, AJ, 
113, 1652 

Gebhardt, K., Kissler-Patig, M. 1999, AJ, October 
issue 

Geisler, D., Lee, M. G., Kim, E. 1996, AJ, 111, 1529 

Gnedin, O.Y., Ostriker, J.P. 1997, ApJ, 474, 223 

Gould, A., Flynn, C, Bahcall, J.N. 1998, ApJ, 503, 
798 

Gratton, R.G., Fusi Pecci, F., Carretta, E., Clcmen- 
tini, G., Corsi, C.E., Lattanzi, M. 1997, ApJ, 491, 
749 

Harris, G.H.L., Harris, W.E., Poole, G.B. 1999, AJ 
117, 855 

Harris, W.E. 1996, AJ 112, 1487, 



Kissler-Patig, M., Gebhardt, K. 1998, AJ, 116, 2237 
Krist, J., Hook, R. 1997, TinyTim v4.4 handbook, 



ittp : //physun. physics .mcmaster . ca/Globular .html 



Harris, W.E. 1997, The Extragalactic Distance Scale, 
Symp.Ser.10 STScI, Cambridge University Press 

Harris, W.E, van den Bergh, S. 1981, AJ, 86, 1627 

Harris, W.E., Allwright, J.W.B., Pritchet, C.J., van 
den Bergh, S. 1991, ApJS, 76, 115 

Harris, W.E., Harris, G.L.H., McLaughlin, D.E. 1998, 
AJ, 115, 1801 

Hilker, M., Kissler-Patig, M., Richtler, T., Infante, 
L., Quintana, H. 1999a, A&AS, 134, 59 

Hilker, M., Infante, L., Vieira, C, Kissler-Patig, M., 
Richtler, T. 1999b, A&AS, 134, 75 

Holtzman J., et al. 1995, PASP 107, 1065 

Jacoby, G.H., Ciardullo, R., Ford, H.C. 1990, ApJ, 
356, 332 

Kauffamnn, G. 1996, MNRAS 281, 487 

Kavelaars, J.J., Harris, W.E., Hanes, D.A., Hesser, 
J.E., Pritchet, C.J. 1999, AJ in press 

Kissler-Patig, M., et al. 1998a, AJ, 115, 105 

Kissler-Patig, M., Richtler, T., Storm, J., Delia Valle, 
M. 1997, A&A, 327, 503 

Kissler-Patig, M., Forbes, D.A., Minniti, D. 1998b, 
MNRAS, 298, 1123 



http : //scivax . stsci . edu/~krist/t inytim.html 



Kundu, A., Whitmore, B.C. 1998, AJ, 116, 2841 

Kundu, A., Whitmore, B.C., Sparks, W.B., Mac- 
chetto, F.D., Zepf, S.E., Ashman, K.M. 1999, ApJ, 
513, 733 

Kurth, O.M., Fritze-v. Alvensleben, U. 1999, A&A 
in press, 



astro-ph 9803192 



Larsen, S.S., Richtler, T. 1999, A&A, 345, 59 

Laurent-Muehleisen S.A., Kollgaard R.I., Ryan P. J., 
Feigelson E.D., Brinkmann W., Siebert J. 1997, 
A&AS, 122, 235 

Lee, M.G., Geisler, D. 1993, AJ, 106, 493 

Lee, M.G., Kim, E., Geisler, D. 1998, AJ, 115, 947 

Maraston, C. 1998, MNRAS, 300, 872 



McLaughlin, D.E. 1999, AJ, in press, |astro-ph/9901283 

McLaughlin, D.E., Pudritz, R.E. 1996, ApJ, 457, 578 

Miller, G.E., Scalo, J.M. 1979, ApJS, 41, 513 

Neilsen Jr., E.H. 1999, PhD thesis, John Hopkins Uni- 
versity 

Ortolani, S., et al. 1995, Nature, 377, 701 

Pierce, M., et al. 1994, Nature, 371, 385 

Poulain, P. 1988, A&AS, 72, 215 

Rcnzini, A. 1999, STScI workshop, When and How do 
Bulges form?, eds. CM. Carollo, H.C. Ferguson, 
R.F.G Wyse, Cambridge University Press 

Rosenberg, A., Saviane, I., Piotto, C, Aparicio, A. 
1999, AJ, November issue 

Saha, A., et al. 1996, ApJ, 466, 55 

Salpeter, E.E. 1955, ApJ, 121, 161 

Scalo, J.M. 1986, Fundamentals of Cosmic Physics, 
Vol. 11, 1 

Schlegel, D.J., Finkbeiner, D.P., Davis, M. 1998, ApJ, 
500, 525 

Schweizer, F. 1997, ASP Conf. Ser. Vol. 116, The 
nature of elliptical galaxies, eds. M. Arnaboldi, 
G.S. Da Costa, P. Saha 

Schweizer, F., Miller, B.W., Whitmore, B.C., Fall, 
S.M. 1996, AJ, 112, 1839 

Seeker, J. 1992, AJ, 104, 1472 



17 



Sharpies, R.M., Zcpf, S.E., Bridges, T.J., Hanes, 
D.A., Carter, D., Ashman, K.M., Geisler, D. 1998, 
AJ, 115, 2337 

Silverman, B.W. 1986, Density Estimation for Statis- 
tics and Data Analysis, Monographs on Statistics 
and Applied Probability, Vol. 26, Chapman & Hall 

Tonry, J.L., Ajhar, E.A., Luppino, G.A. 1990, AJ, 
100, 1416 

van den Bergh, S. 1994, AJ, 108, 2145 

de Vaucouleurs C, et al. 1991, Third Reference Cat- 
alogue of Bright Galaxies, Version 3.9 

Vesperini, E., Heggic, D.C. 1997, MNRAS, 289, 898 

Whitmore, B.C. 1997, The Extragalactic Distance 
Scale, Symp.Ser.10 STScI, Cambridge University 
Press 

Williams, R.E., et al. 1996, AJ, 112, 1335 

Worthcy, G. 1994, ApJS, 95, 107 

Zepf, S.E., Ashman, K.M. 1993, MNRAS, 264, 611 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v4.0. 



18 



Table 1 

General properties of NGC 4472 



NGC 4472 



RA(2000) a 12h 29m 46.79s ±1.25s 

DEC(2000) a +08° 00' 01.50" ±1.25' 

m v h 8.41 ± 0.06 mag 

V-I c 1.24 ±0.01 mag 

B-V b 0.96 ±0.01 mag 

MType b E2/S0(2) 

Vhdio b 868 ± 8 km s" 1 

Distance d 15.8 ± 0.8 Mpc 



a: Laurent-Muchleiscn S.A., et al. 1997 
b: RC3v9 (de Vaucoulcurs, et al. 1991) 
c: Poulain P., 1988 
d: this paper 



Table 2 
Summary of the observations 



Prg. ID + PI 


RA(2000) 


DEC(2000) 


F555W exp. time a 


F814W exp. time a 


GO. 5236 Westphal 


12 29 48.73 


+07 59 55.85 


1800 sec 


1800 sec 


GO. 5920 Brodic (north) 


12 29 45.22 


+08 02 34.38 


2200 sec 


2300 sec 


GO. 5920 Brodic (south) 


12 29 45.20 


+07 57 30.11 


2200 sec 


2300 sec 



a Effective exposure time, after images were stacked. 
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Table 3 

Turn-over magnitudes and dispersions for each sub-population luminosity function 



Filter 


Population 


M TO 


AM TO 


°"i 5 






Total sample 














V 


blue 


23.62 


0.09 










red 


24.13 


0.11 








I 


blue 


22.48 


0.07 










red 


22.90 


0.11 








Westphal (inner) sample 














V 


blue 


23.55 


0.15 


1.21 


0.10 


118 




red 


23.97 


0.19 


1.42 


0.10 


180 


I 


blue 


22.62 


0.21 


1.24 


0.14 


129 




red 


22.83 


0.19 


1.39 


0.10 


174 


Brodie (outer) sample 














V 


blue 


23.68 


0.11 


1.36 


0.09 


175 




red 


24.30 


0.13 


1.39 


0.09 


232 


I 


blue 


22.33 


0.08 


1.14 


0.07 


178 




red 


22.97 


0.13 


1.39 


0.08 


213 



The TO magnitudes are neither corrected for reddening, nor for the aperture correction discussed in Sect. 2.2.3 



The corrections for galactic reddening are Ay = 0.069 and Ai = 0.040. The corrections for aperture, in addition 
to the correction of Holtzman C = 0.1 (which is already included) are Cy = 0.050 ± 0.015 and Ci = 0.080 ± 0.010 
(which are not included here). 
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Aperture diameter [pix] 



Fig. 1. — Curves of growth for PSF models and selected data. The curves represent different models (solid line 
= raw TinyTim PSF; long-dashed line = Hubble law (0.05" core radius) ® PSF; short-dashed line = Hubble law 
(0.1" core radius) ® PSF; dot-dashed line = Hubble law (0.2" core radius) ® PSF while the open circles show the 
selected globular cluster data: 13 candidates in V and 15 in I, 
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Fig. 2. — Completeness curves for the central pointing (dashed lines) and for the north- and south-pointings (solid 
lines). The central images have much higher background values than the northern and southern pointings. Each 
curve refers to a particular background value bin, i.e. from faint to bright magnitudes in counts: 0-50, 50-100, 
150-200, 200-250, 250-350, 350-550. 
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Fig. 3. — Background contamination histograms for the Vf606VK and Ifshw band resulting from the HDF-N. The 
unshaded histograms show the background object luminosity distribution as obtained from HDF-N directly, while 
the shaded ones show the same data after applying our criteria. The lines are exponential least square fits to the 
small histograms. In order to ensure a reasonable fit we set the upper limit for the magnitude range at the 50% 
completeness level quantified during the completeness tests. Using the exponential law, the integrated number of 
background objects is ~ 30 up to the 50% completeness limit in both filters. 
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. 4. — Color magnitude diagram of globular clusters in NGC 4472. The box delineates the initial cuts made to 
raw data prior analysis. 705 globular cluster candidates are enclosed. 
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Fig. 5. — Color distribution of globular clusters in NGC 4472. A bimodal distribution is found to be present at 
the 99.99% confidence level using the KMM mixture-modeling algorithm (Ashman et al. 1994). The two peaks 
lie at 0.99 ± 0.01 and 1.24 ± 0.01 (0.93 ± 0.01 and 1.18 ± 0.01 corrected for reddening and aperture). A further 
analysis using a density estimation technique (kernel estimator with an Epanechnikov kernel) yields color peaks at 
V — I = 1.00 ± 0.01 and 1.22 ± 0.01. The solid line depicts the derived distribution function with statistical upper 
and lower error limits (short-dashed lines). 
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Fig. 6. — Globular cluster luminosity functions of the halo (Brodie) pointings (upper 4 plots) and central (Westphal) 
pointing (lower 4 plots). The number of clusters that went into the histograms up to the faintest completeness 
limit are indicated. The histograms are shown for the selected globular clusters, uncorrected for completeness and 
background variations (see text). The dashed curves are the best fitting t$ functions as derived from our maximum 
likelihood analysis, taking into account completeness and background corrections. The t$ functions were arbitrarily 
normalized. The vertical dotted lines mark our 50% completeness limit at the lowest background value. 
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Fig. 7. — Radial distributions of the red (open circles) and blue (stars) globular clusters. The dashed line represents 
the best weighted linear least-square fit to the red population. The dot-dashed line gives the best fit to the blue 
population. The solid line shows the best fit for the entire globular cluster population. The gradient is twice as 
large as for each sub-population. 
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Fig. 8. — Upper panel: Radial globular cluster size distribution for blue (open stars) and red (open circles) globular 
clusters. The dot-dashed line is the best least-square fit to the blue sub-population while the dashed line gives 
the best fit to the red sub-population. Lower panel: Histogram of relative globular cluster sizes for the red 
population (solid histogram) and the blue population (dashed histogram). The sizes of the blue globular clusters 
are substantially larger than those of the red. 
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Fig. 9. — Maraston's (1998) Simple Stellar Population models. Upper panels: I-band data. Lower panels: V-band 
data. The solid grid lines represent iso-metallicity tracks of Z — 0.001, 0.006, 0.02, 0.04. Dash-dotted lines are 
isochrones of 9, 10, 11, 13, 14 and 15 Gyr. Solid error bars represent la errors and dashed error bars refer to 3cr 
errors. The averaged maximum age difference (age D i U c— age ro d) is AT = 0.7 ± 1.8 Gyr. 
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Fig. 10. — Simple Stellar Population models by Kurth et al. (1999) calculated for a single stellar burst of 10 7 
Gyr duration. Upper panels: I-band TO magnitudes. Lower panels: V-band TO magnitudes. Within the grids 
the solid lines (vertical) represent iso-metallicity tracks of Z = 0.0001, 0.0004, 0.004, 0.008, 0.02, 0.05 (from left 
to right) which correspond to [Fe/H]= —2.3, —1.7, —0.7, —0.4, 0.0 and 0.4. The dash-dotted (horizontal) lines are 
isochrones of 8, 10, 12, 14 and 16 Gyr. The solid error bars represent the la errors. The dotted error bars give the 
3<7 deviations. An averaged maximum age difference of AT = — 2.3 ± 1.6 Gyr is derived. 
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Fig. 11. — Simple Stellar Population models by Bruzual & Chariot (1996). Upper panels: I-band data. Lower 
panels: V-band data. The solid grid lines represent iso-metallicity tracks of Z = 0.0004, 0.004, 0.008, 0.02 and 
0.05. Dash-dotted lines are isochrones of 8, 10, 12, 14, 16 Gyr. Using these grids an averaged maximum-age 
difference of AT = 3.3 ± 2.2 Gyr was derived. 
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Fig. 12. — Worthey's (1994) Simple Stellar Population models. Upper panels: I-band data. Lower panels: V-band 
data. The solid grid lines represent iso-metallicity tracks of Z = 0.0002, 0.0004, 0.004, 0.008, and 0.02. Dash- 
dotted lines are isochrones of 8, 10, 12, 14, 16 Gyr. The averaged maximum age difference is AT = —4.0 ± 1.6 
Gyr. 
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